Cavity sideband cooling of trapped molecules 
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The efficiency of cavity sideband cooling of trapped molecules is theoretically investigated for the 
case where the IR transition between two rovibrational states is used as a cycling transition. The 
molecules are assumed to be trapped either by a radio-frequency or optical trapping potential, de- 
pending on whether they are charged or neutral, and confined inside a high-finesse optical resonator 
which enhances radiative emission into the cavity mode. Using realistic experimental parameters 
and COS as a representative molecular example, we show that in this setup cooling to the trap 
ground state is feasible. 

PACS numbers: 37.10.Mn,37.30.+i 



I. INTRODUCTION 



With the field of cold and ultracold molecules coming 
of age in the past few years, efficient cooling methods 
are becoming more and more valuable. Optical methods 
have the advantage that they do not rely on collisional 
processes which are hard to predict and which often open 
severe collisional loss channels. A number of ingenious 
optical cooling schemes for molecules have been devised: 
It has been realized, e.g., that translational cooling of 
molecules might be possible with light sources with a 
specifically designed set of lines [l], Q . Alternatively, the 
almost-closed level scheme a molecule like CaF has been 
exploited @ to laser cool molecules as it is usually done 
with Q • Even a 3D scheme has been proposed resembling 
the popular magneto-optical trap for atoms Q . Zcppen- 
feld et al. have proposed to use electro-optical forces to 
carry away up to IK per spontaneously decayinginfrared 
photon [f|. Optical cavities have been predicted @,[3 and 
used for atoms [§] as a way to achieve translational cool- 
ing without relying on spontaneous emission , in principle 
allowing cooling of molecules. In Refs. (Tfj|.]Tl| we have 
shown that cavity-enhanced Raman scattering can enable 
one to cool both the internal and the external degrees 
of freedom of molecules, using simulations based on ab- 
initio calculations for OH and NO radicals. In this pro- 
posal the cavity enhances the anti-Stokes, cooling tran- 
sitions over other heating transitions. In a recent paper 
Vuletic and co workers demonstrated sideband cooling of 
the motion of a single trapped ion to the ground state 
of the trap potential via cavity enhanced scattering fl2| . 
In the present article we make the next step and demon- 
strate theoretically cavity-enhanced sideband cooling of 
molecules and molecular ions in a strongly confining ex- 
ternal potential, in the setup sketched in Fig. [1] . We 
show that cooling to the ground state of the external po- 
tential can be achieved for experimentally feasible setups. 
The new scheme thus offers the possibility to reach lower 
final temperatures compared to the previously reported 



cooling schemes 0, [T(| , where the limit is given by the 
cavity linewidth. 

This article is organized as follows. In Sec. [II] the theo- 
retical model is introduced and the derivation of the rate 
equations for cooling the external motion of the molecules 
is sketched. Moreover, the method used for obtaining the 
numerical results is described. In Sec. IIIII the results are 
reported and discussed, while the conclusions are drawn 
in Sec. HV1 



II. THE MODEL 

We theoretically explore the efficiency of performing 
sideband cooling of a molecule of mass M, which is con- 
fined by a harmonic potential and couples to the mode 
of a high-finesse resonator. The level scheme is sketched 
in Fig. [21 a two-level transition of the molecule is driven 
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FIG. 1: Schematic representation of the setup. A molecule 
is confined by an external harmonic potential inside a high- 
finesse resonator. The transition between two molecular vi- 
brational levels is driven transversally by a laser and scat- 
ters photons into a quasi-resonant cavity mode. The scatter- 
ing processes are tailored in order to enhance the transitions 
which cool the external motion. 
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transversally by a laser with coupling strength il (Rabi 
frequency), and emits photons at frequency u> c . The bad- 
cavity limit is assumed, which implies that the molecule 
scatters prevalently photons into the cavity mode, which 
are then lost via cavity decay. 

Let us now provide a specific example. The two-level 
transition, which is driven by electromagnetic radiation, 
is composed of two vibrational states v = 0, v = 1 of the 
electronic ground state of a molecule, which we denote 
by \g) and |e), respectively. Let coo be the transition fre- 
quency. This level scheme is closed inside the vibrational 
subspace. Spontaneous emission into different rotational 
substates is suppressed by the cavity, while any residual 
leakage can be handled by microwave re-pump schemes. 
These properties allow one to use this transition for cav- 
ity cooling. 

We further note that this sort of vibrational transitions 
is usually in the infrared regime (100-4000 cm -1 ) and has 
typically a homogeneous lincwidth below 1 kHz, which 
allows one to spectrally resolve the motional sidebands. 
In order to obtain convenient cooling rates, and hence 
efficient cooling, spontaneous emission is here enhanced 
by coupling with the mode of an infrared high-finesse 
resonator. The spatial confinement can be obtained by 
means of dipole or ion traps, depending on whether the 
molecule is neutral or charged fl3Ml5l ] . 

In order to perform sideband cooling, the initial spa- 
tial localization of the molecule must be smaller than the 
wavelength of the driving laser. This requirement can 
be rewritten as r/y2(n) + 1 <C 1, with (n) the mean oc- 
cupation of phononic states of the harmonic oscillator 
and rj the Lamb-Dickc parameter, r) = ^Jhk 2 /2Mv, with 
k the wave vector of the transition and v the trap fre- 
quency in angular units[20l.[2ll]. One necessary condition 
is that the Lamb-Dicke parameter is smaller than unity 
i]€l. Moreover, the molecule should be pre-cooled, e.g., 
by means of an optical cooling schemes [3,111 or other non- 
optical methods @. Nevertheless we note that sideband 
cooling can be efficient even for r\ ~ 1 Moreover, 
schemes for ground-state cooling for Lamb-Dicke param- 
eters larger than unity have been theoretically proposed 
in fill ] , and can be extended to situations in which dissi- 
pation occurs via cavity decay. 



A. Master equation for cavity sideband cooling 

We now introduce the model, from which we extract 
the cooling rates and the final temperature expected by 
cavity sideband cooling the molecules. Below we assume 
the Lamb-Dicke regime, and hence expand in first order 
in the small parameter rj (2^, [2l| . The model has been 
extensively reported, for instance in Refs. fl7l. Il8j. 

We study the coupled dynamics of the center-of-mass 
motion of the molecule, the two-level transition, and the 
cavity mode in the reference frame rotating at the fre- 
quency lul of the laser, driving the two-level molecular 
transition. Denoting by p the density matrix for molcc- 
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FIG. 2: Schematic representation of the scattering process 
leading to cooling. The molecular vibrational levels are la- 
beled by v — 0, 1, the phononic, trap excitations by n. The 
pump laser at frequency ujl drives resonantly the red side- 
band transition \v = 0, n) — » \v = l,n — 1). The laser photon 
is scattered into a cavity mode at frequency uj c = ujq, which 
couples resonantly with the molecular transition. Since the 
photon is then lost via cavity decay, the molecule is cooled by 
one phonon. 

ular internal and external degrees of freedom (in this re- 
duced Hilbert space) and for the cavity mode, its dynam- 
ics is governed by the master equation 

d t p= -^[H,p} + C K p + C y p, (1) 

where the Hamiltonian H accounts for the coherent evo- 
lution, and the supcropcrators C K and £ 7 describe cavity 
decay and spontaneous emission, respectively. We now 
give the specific form of each term. The full Hamiltonian 
can be written as the sum of the terms 

H = H M +H T + H C + V I , (2) 

where 

H M = -HA\e){e\ (3) 

is the Hamiltonian for the internal degrees of freedom, 
with A = uj^—loq. The Hamiltonian Ht describes the dy- 
namics of the harmonic oscillator for the center-of-mass 
motion, which is here assumed to be one dimensional for 
convenience, 

H T = £ + ±Mv>X' = *,(lfib + ±), (4) 

with X and P position and canonical conjugate momen- 
tum of the molecule, and b = (iP + MvX) / y/2hmv an- 
nihilates a quantum of motional excitations hv such that 
[6,6'''] = 1. We denote by {\n)} the number states with 
n = 0,1,2,..., which are eigenstates of the Hamiltonian 
in Eq. Q. The Hamiltonian for the cavity mode is 

He = —hS c a'a , (5) 

where operators a and annihilate and create, respec- 
tively, a cavity photon at frequency oj c and 5 C = lu^ — uj c 
is the detuning between laser and cavity mode. 
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In the bad-cavity regime we may reduce the Hilbcrt 
space of the cavity mode to the states |0 C ) and |l c ), 
namely, zero and one photon in the cavity mode, respec- 
tively. The coupling of the molecule with the laser and 
with the cavity mode are given by operator Vi which 
is here reported in the Rotating Wave Approximation 
(RWA), as the pump is assumed to be sufficiently weak 
to neglect the effect of counter-rotating terms. The in- 
teractions in the Lamb-Dickc expansion are given by 



Vi = Vi 



(o) 



V^ib + b^ + Oij] 2 ), 



where the superscript gives the correspon ding perturba- 
tive order. The zero-order term reads fl7l Il8j ] 



Vj {0) = h(il + go) cos <i>)o + H.c. 



(6) 



where fl is the Rabi frequency, g the vacuum Rabi fre- 
quency and 4> the phase between a maximum of the cav- 
ity standing-wave and the center of the trap. Operators 
a = \g)(e\ and its adjoint crt flip the molecular excitation. 

The term in first order in the Lamb-Dicke expansion 
include the coupling with the external degrees of freedom 
and read [U EH 

Vj 1 ' 1 = hr](iQ cos &l — go* cos C sin <j))a + H.c. , (7) 

where Ol (&c) is the angle between the axis of the mo- 
tion and the laser (cavity) wave vectors. 

The superoperators, describing dissipation in Eq. (fTJ, 
take the form 



£ K P = -x {Zapa 1 - {a^a,p}) , 



£~«P = o (2crpcr f - {crV, p}) 



(8) 
(9) 



with the molecular (cavity) linewidth 7 (re) and p ac- 
counting for the mechanical effects of the spontaneously 
emitted photon on the motion (l7j . Since spontaneous 
decay is very slow with respect to the other rates, in the 
simulations we will take p ~ p, hence discarding terms 
of the order r] 2 j which describe the diffusion due to the 
recoil of the spontaneously emitted photons. 



B. Rate equations 

We characterize the efficiency of cavity-enhanced side- 
band cooling by means of a rate equation, which gives the 
dynamics of the occupation of the oscillator levels \n), 
p n = Tr{\n)(n\p}. This treatment is valid in the limit 
in which the dynamics of internal and cavity degrees of 
freedom is faster than the external motion, which requires 
i|< 1. In this case, one can assume that the coupled sys- 
tem, given by the cavity mode and two-level transition, 
reaches a steady state given by the density matrix g$, 
which is the solution of the equation CqQs = 0. Here, 
Qs is defined over the Hilbert space of cavity mode and 
internal degrees of freedom of the molecule and 
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[H M + H c + V I (Q \g}+C^g + C K g (10) 



is the Liouvillc operator for internal and cavity degree 
of freedom. Under these assumptions the rate equation 
determining the dynamics of the trap states takes the 
form 

i>n = -{nA- + (n+l)A + )p n + (n+l)A-p n+ i + nA + p n -i , 
with fij 



A± 



-2Re 



Triy« 



1 



C T w 



V^gs 



(12) 



Here, A + is the heating (cooling) rate and Vj 

is defined in Eq. ([?])• These are found by assuming 
that the Rabi frequency Q is the smallest parameter, 
Q <C g,j,K,u, which implies that the relevant internal 
and cavity states which are involved in the dynamics 
are the ground state \v = 0,0 C ) and the excited states 



0,l c 



1,0 C ). In the bad-cavity limit, state 



\v = 0, l c ) decays very fast, while \v = 1,0 C ) is charac- 
terized by a relatively smaller rate of decay, such that 
k g >• 7 pjl EH ■ An effective cooling rate W can be 
defined for the case A- > A + [20l |. 



W = A- 



(13) 



which gives the rate at which phonons are exponentially 
damped towards the steady-state value of the trap. The 
steady-state expectation value of the phonon occupation 
number is given by [20| 



(n)st 



A, 



A- - A 4 



(14) 



C. Numerical evaluation of the master equation 

We check the analytical predictions by means of nu- 
merical calculations, which consists in performing the 
time evolution of the density matrix of cavity and molec- 
ule internal and external degrees of freedom according to 
Eq. (fTJ . This is performed taking the Hilbcrt space of the 
cavity mode to be reduced to the photon states |0 C ), |l c ), 
which is a reasonable choice considering that in the bad- 
cavity limit the cavity mode is either in the vacuum or 
has an occupation of one photon. Moreover, the molec- 
ular vibrational levels are approximated by a two-level 
system, which is resonantly driven by radiation. Finally, 
the center-of-mass oscillator is approximated by the first 
5 levels. A convergence check for an increasing number 
of trap states has been carried out to ensure that the 
error due to the truncation does not significantly affect 
the results. The resulting density matrix has dimensions 
20x20. The integration of Eq. (fTJ) is done with an fourth- 
order Runge-Kutta scheme to obtain the time evolution 
of the coupled quantum system. The time evolution of 
the center-of-mass motion is found by taking the partial 
trace of the total density matrix, as a function of time, 
over the internal and cavity degrees of freedom. In or- 
der to determine the cooling rate, it is assumed that the 
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TABLE I: Parameters of the chosen molecules and their IR 
cooling transitions. The point group of the molecule and 
the irreducible representation of the corresponding vibrational 
normal mode is given. The transition frequencies, transition 
dipole moments(lau = 2.54Debye) and Einstein A coeffi- 
cients 7 are approximated from the ab-initio frequency anal- 
ysis. 

PG Irred. vn [cm -1 ] /^i2[au] 7[s _1 ] 



CHBr 3 

HCCCF 3 

TMA 

COS 

CFI 3 

CSCla 

MgH + 



Gzv E 
Czv A\ 
Czv Ai 

Coov ^g 

Czv Ai 

C2v A\ 



635 
1234 
2910 
2108 
1038 
1131 
1609 



0.10 

0.15 
0.067 

0.15 
0.087 

0.13 
0.036 



5.(i 
80 
226 
424 
17.1 
5.12 
11 



rate equations, Eq. (jlip . apply. Thus the time evolution 
of the trap states is fitted against the populations p n to 
obtain A + and A_, which allows us then to extract W 
and (n) st- 



ill. RESULTS 

In order to evaluate the efficiency of cavity sideband 
cooling of molecules, we first identify some good candi- 
dates. We note that a strong vibrational transition favors 
large scattering into the cavity. Moreover, for the case 
of neutral molecules trapped by a far-detuned standing- 
wave laser a large polarizability is needed to generate a 
sufficiently deep potential well. Several molecular can- 
didates have been selected and calculated with the soft- 
ware package Gaussian (22j . The equilibrium geometry of 
the molecules has been optimized with density functional 
theory and a hig h-quality atomic basis set (B3LYP/aug- 
cc-pVTZ [23l - |25| L Based on these geometries a frequency 
analysis and polarizability calculation have been carried 
out to determine the spontaneous emission rates 7 of the 
excited vibrational state and polarizability tensors as the 
decisive molecular properties. A small choice of possible 
candidates which we have considered for implementing 
cavity sideband cooling is shown in Tab. |IJ 

The molecule COS is selected from the list as a neu- 
tral candidate for its remarkable strong IR transition in 
the asymmetric stretch vibration. Its large polarizabil- 
ity makes it also a good choice for optical trapping. For 
the two-level system v — and v = 1 of the asymmet- 
ric stretch mode with the parameters given in Tab. U arc 
chosen. The test candidate COS is assumed to reside in 
a standing- wave optical potential which has a trap depth 
of ss 900 \iK and a trap frequency v = 2ir x 350 kHz 
in the harmonic approximation. Such a potential could 
for example be formed by a laser with a wavelength of 
532 nm and an effective power of 500 W achieved by a 
build-up cavity. The resulting Lamb-Dicke parameter is 
then 77 = 0.02. Molecular samples below lmK may be 



achieved in the near future [6j, making it plausible to 
assume the Lamb-Dicke regime in our system. 

The cavity axis and the pump laser axis is chosen to 
be in a 45° arrangement with the trap axis (©c = ©l = 
7r/4). The amplitude of the electric field of the cavity 
mode is e c = 150 Vm -1 and the corresponding cavity 
line width is k = 2tt x 5 MHz. The single-particle coop- 
erativity C\ = <? 2 / K 7 = 61 warrants the strong-coupling 
regime which is needed to provide a sufficiently strong 
enhancement of the scattering rate. However, it is im- 
portant to notice that the cavity coupling g is smaller 
than the trap frequency (g = 0.41 v). This is a major 
difference to previous theoretical works [TH, [l|| which 
assumed significantly larger values of g with respect to 
i>, and which determines in our case the di sapp earance 
of most interference phenomena reported in [l8[ 0j| . 

In order to explore the parameter regime identified for 
the case of the COS molecule, the cooling and heating 
rates are calculated by the two different methods de- 
scribed in Sec. Ill CI namely, the analytical predictions 
from perturbation theory and the numerical results found 
by integrating master equation Eq. ([I]). 

Figure [3] displays the contour plots of cooling rate 
and steady state phonon number as a function of laser- 
molecule detuning A and the laser-cavity detuning 6 C . 
The top row in Fig. [3] shows to the predictions of the 
analytical model, while the second row reports the corre- 
sponding predictions of the numerical results for the same 
set of parameters. Their comparison shows discrepancies 
between analytics and numerics. This is understood: In 
fact, the pcrturbative treatment is based on the assump- 
tion that the Rabi frequency 17 is the smallest parame- 
ter, and that the excited state of the two-level transition 
is essentially empty [l8j, while for the parameters here 
considered with 7 being the smallest rate this is not ful- 
filled. We have also checked that in the limit where the 
perturbative treatment is expected to be valid, the pre- 
dictions of analytics and numerics are almost identical. 
In the following we refer then to the numerical results. 
The analytical results, nevertheless, serve for a double 
purpose. In first place, they allow one to identify the 
parameter interval where cooling is optimal (allowing to 
quickly find the region where the numerical effort should 
be concentrated on). Moreover, since the dynamics de- 
scribed by the analytical solution is well understood, see 
e-g- 13, 13, by comparison with the numerical results 
they permit one to identify the physical origin of the dis- 
crepancies. 

We first notice that efficient cooling is found for the 
values of the detuning at which the red sideband of the 
narrow dressed state of cavity and molecule is driven. At 
large values of the detuning |<5 C | — > 00, when the coupling 
with the cavity mode becomes irrelevant, this dressed 
state tends to the molecular vibrational state v = 1 and 
optimal cooling is thus found when A = —v, which is the 
sideband cooling condition in free space [2(], [2l| . Cor- 
respondingly, the temperature is minimal. We also find 
that the numerical results predict that the cooling rate 
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FIG. 3: Contour plots of the cooling rate (left column) and corresponding steady-state occupation (right column) of the 
motion of a trapped COS molecule driven at the asymmetric stretch vibration (2108 cm" 1 ) as a function of the laser- molecule 
detuning A and the laser-cavity detuning 8 C . The parameters are: g = 0.41^, r\ = 0.02, k = lAu, 7 = 1.9 ■ 10~ 4 v and 
Ol = 9c = <t> = 7r/4. Subplots a) and b) correspond to the results found by applying perturbation theory Eq. (|12[) and taking 
Q, — 0.05 v. Subplots (c) and (d) are found for the same parameters as in (a) and (b), but are calculated from the density 
matrix found by numerical integration of Eq. (jXJ . Subplots (e) and (f) are the same as (c) and (d) but with pump strength 

n = 0.1 ia 



exhibits a maximum (and the finite occupation a mini- 
mum) about a finite value of |5 C |. Comparison with the 
analytical results shows that this behavior must be due 
to saturation effects, which are not accounted for in the 
perturbative treatment. Indeed, for the chosen parame- 
ters the laser saturates the red sideband transition (be- 
ing vfcl > 7). As a consequence, the mean vibrational 
occupation at steady state, extracted from the numerical 
method and displayed in Fig. |3|d), shows oscillations as 
the parameters are varied across the sideband resonance. 



The cooling rate which we extract at the maximum is 
W = 2 • 10 3 s _1 . 

Figure |3|e) and Fig.[3ff) display the numerical predic- 
tions for cooling rate and steady state occupation when 
the parameters are the same as in Figs. [3{c) and Fig. 
|3{d), but the strength of the pump is doubled, = 0.1 v. 
Larger cooling rates are achieved in this case, and the 
minimal temperature is also increased. The observable 
power broadening leads also to a larger parameter space 
where efficient cooling is observed. 
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Q. [units of v] 

FIG. 4: Extremal values of the cooling rate (black line) and of 
the steady state occupation (grey line) versus the laser Rabi 
frequency (in units of v). The other parameters are the same 
as in Fig. [3] 

Figure^] shows the maximum cooling rate as a func- 
tion of VL, beyond the validity of the perturbative regime, 
which predicts that the cooling rate scales quadratically 
with [lH. We observe that with increasing pump 
strengths the cooling rate becomes larger, and increases 
almost linearly with Q. Correspondingly, the steady- 
state phonon number increases exponentially with f2. 
Thus the choice of the pump strength is a trade-off be- 
tween final temperature and fast cooling process. 

IV. CONCLUSION 

In this article we have given numerical evidence that 
cavity-enhanced sideband cooling of trapped molecules is 



feasible, by performing the calculations for a particular 
molecule, COS. The chosen example of COS has favor- 
able properties such as a large polarizability and a strong 
vibrational transition. We have identified experimental 
parameters for trapping and cooling the COS molecule, 
which may be realized in the near future. The use of 
infrared transitions well-resolved sidebands may easily 
be achieved in different types of harmonic traps. As 
long as the strong-coupling regime can be implemented 
with a chosen vibrationally transition, concurrent decay 
channels are outrun by the enhancement. Moreover the 
scheme is well suited for ion traps due to their long stor- 
age times and their low anharmonicities. Under these 
considerations ground-state cooling on a sub-millisecond 
time scale seems possible. More efficient cooling could be 
achieved by increasing the the cavity coupling g, which 
could be achieved by exploiting collective enhancement 
in a molecular gas as proposed in Ref. (26[. Further im- 
provement of the cooling process is conceivable by adap- 
tive control of the Rabi frequency in time. 
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